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Abstract 

We investigate superconductor/insulator/ferromagnet/superconductor (SIFS) tunnel Josephson junctions in the dirty limit, using 
the quasiclassical theory. We consider the case of a strong tunnel barrier such that the left S layer and the right FS bilayer are 
decoupled. We calculate quantitatively the density of states (DOS) in the FS bilayer for arbitrary length of the ferromagnetic 
layer, using a self-consistent numerical method. We compare these results with a known analytical DOS approximation, which 
is valid when the ferromagnetic layer is long enough. Finally we calculate quantitatively the current-voltage characteristics of a 
SIFS junction. 
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1. Introduction 

It is well known that superconductivity and ferro- 
magnetism are two competing orders. The coexistence 
of singlet superconductivity and ferromagnetism is ba- 
sically impossible in the same compound but may be 
easily achieved in artificially fabricated superconduc- 
tor/ferromagnet (S/F) hybrid structures. In this case, the 
coexistence of the two orders is due to the proximity 
effect [1,2]. The main manifestation of the proximity 
effect in S/F structures is the damped oscillatory behav- 
ior of superconducting correlations in the F layer. Two 
characteristic lengths of the decay and the oscillations 
are, respectively, ^/i and E,f2- 

Therefore in S/F heterostructures there is a unique 
possibility to study the properties of superconducting 
electrons under the influence of the exchange field in 
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the ferromagnet. Recent progress in the preparation of 
the high quality S/F layered structures permitted to 
experimentally observe many striking phenomena that 
are quite interesting for applications such as a non- 
monotonic dependence of their critical temperature and 
oscillations of critical current in S/F/S Josephson junc- 
tions as a function of the F layer thickness (see [1] and 
references therein). It is possible to fabricate Joseph- 
son % junctions with a 7T-phase difference in the ground 
state, which are good candidates for elements in super- 
conducting logic circuits [3]. 

SIFS junctions, i.e. S/F/S trilayers with one transpar- 
ent interface and one tunnel barrier between S and F 
layers, represent an interesting case of K junctions for 
applications where active Josephson junctions are re- 
quired. The SIFS structure offers the freedom to tune 
the critical current density over a wide range and at the 
same time to realize high values of the product of the 
junction critical current I c and its normal state resis- 
tance Rn [4]. In addition, Nb based tunnel junctions are 
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Fig. 1. Geometry of the considered system. The thickness of the 
ferromagnetic interlay er is df. The transparency of the left S/F 
interface is characterized by the coefficient Jg\ , and the transparency 
of the right F/S interface is characterized by y S 2- 

usually underdamped, which is desired for many ap- 
plications. SIFS junctions are also interesting from the 
fundamental point of view since they provide a conve- 
nient model system for a comparative study between 
0-7T transitions observed from the critical current and 
from the density of states (DOS) [5]. 

The purpose of this work is to provide a quantitative 
model describing the DOS in SIFS junctions and to cal- 
culate current- voltage characteristics of a SIFS junction. 
The latter may be used for estimation of the dissipation 
in SIFS-based qubit systems. 

2. Model and basic equations 

The model of an S/F/S junction we are going to study 
is depicted in Fig. 1 and consists of a ferromagnetic 
layer of thickness df and two thick superconducting 
electrodes along the x direction. The left superconduc- 
tor electrode is voltage-biased. Left and right super- 
conductor/ferromagnet interfaces are characterized by 
the dimensionless parameters jg\ and y^, respectively, 
where y B \,B2 = RBi,B20 n /^„, Rb\,b2 are the resistances 
of the left and right S/F interfaces, respectively, a„ is 
the conductivity of the F layer, £,„ — y/Df/2nT c , Df is 
the diffusion coefficient in the ferromagnetic metal and 
T c is the critical temperature of the superconductor (we 
assume h = kg = 1). We also assume that the S/F inter- 
faces are not magnetically active. We will consider the 
diffusive limit, in which the elastic scattering length £ is 
much smaller than the characteristic decay length 

In this paper we concentrate on the case of a SIFS 
tunnel Josephson junction, when Yb\ 3* 1 (tunnel bar- 
rier) and = (fully transparent interface; however in 
our numerical calculations we use finite but very small 
7b2 "C 1). In this case the left S layer and the right FS 
bilayer in Fig. 1 are decoupled. Therefore we can calcu- 
late the quasi-particle current through a SIFS junction 
using the standard formula of the tunnel theory [6], 
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dEN s (E - eV)N{E) [f(E - eV) - f(E)} , 

(1) 



where N S {E) = \E\®(\E\ — A)/VE 2 —A 2 is the BCS 
density of states [&(x) is the Heaviside step function] 
and N(E) is the density of states in the FS bilayer at 
the free boundary of the ferromagnet (x = —df/2). Both 
Ns(E) and N(E) are normalized to their values in the 
normal state; f(E) = [1 + exp(E /T)]~ l is the Fermi 
function. To obtain N(E) we should solve the Usadel 
equations in the ferromagnetic layer of the FS bilayer. 

Using the 0-parameterizations of the normal and 
anomalous Green functions, G = cos 6, F = sin 9, we 
can write the Usadel equations in the F layer as [7,8] 
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where a positive (negative) sign in front of h corresponds 
to the spin up state j" (spin down state \), CO = 2nT(n + 
j) are the Matsubara frequencies, h is the exchange field 
in the ferromagnet, and the parameter x m is the spin- 
flip scattering time. We consider a ferromagnet with 
strong uniaxial anisotropy, in which case the magnetic 
scattering does not couple the spin up and spin down 
electron populations. 

In the S layer the Usadel equations take the form 
(where we omit subscripts 'j" (\)' because equations in 
superconductor look identically for spin up and spin 
down electron states), 

D s d 2 e s 



co sin 9 S — A(x) cos S 



(3) 



2 dx 2 

They should be completed with the self-consistency 
equation, 



A(x)ln- 



v- f2A(x) 
^ \ co sT 
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Here D s is the diffusion coefficient in the superconduc- 
tor and A(x) is the superconducting pair potential. 

Since we calculate the DOS at the free boundary of 
the ferromagnet in the FS bilayer, we need to set to zero 
the 9/ derivative at the left S/F interface, 

(de f /dx)_ df/2 = o. (5) 

At the right F/S interface the boundary conditions are 
given by the relations [9], 



$ n Y(ddf/dx) df/2 = £ v {d9 s /dx) df/2 , 
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where y= ^ V (7„/^„(7 V , a s is the conductivity of the S 
layer and £ s = \/D s /2kT c . The parameter y determines 
the strength of suppression of superconductivity in the 



2 




ElA 2 



ElA 2 




ElA 2 °' 5 L 



ElA 2 




1.05 
















1 
f 


c) 1.00 


J 


1 1 



1.00 



0.95 





1.00 


r 






c) 0.96 



ElA 2 



4 ElA 



ElA 2 




ElA 2 




ElA 8 












1.002 
1.000 




d) 





h) 



ElA 



Fig. 2. (Color online) DOS on the free boundary of the F layer 
in the FS bilayer calculated numerically in the absence of spin-flip 
scattering. Plots (a)-(d) are calculated for h = kT c and (e)-(h) for 
h = 3 %T C . The temperature T = 0.1T C . Parameters of the F/S interface 
are y= y B1 =0.01 (a),(e): d f /$ n =0.5; (b),(f): d f /$„ = 1; (c),(g): 
df/%„ = 2; (d),(h): df/%„ = 3. The approximate analytical solution 
[5] is shown by dashed red lines. 

right S lead near the interface compared to the bulk: no 
suppression occurs for 7=0, while strong suppression 
takes place for y^> 1. In our numerical calculations we 
will assume small y <C 1 . 

To complete the boundary problem we also set a 
boundary condition at x = °o, 

8 s (°o) = arctan=-. (7) 

3. Density of states in the FS bylayer 

To calculate the DOS N(E) on the free F layer bound- 
ary (x — —df/2) we use the self-consistent two-step it- 
erative procedure [10]. In the first step we calculate the 
pair potential's coordinate dependence A(x) using the 
self-consistency equation in the S layer, Eq. (4). Then, 



Fig. 3. (Color online) DOS on the free boundary of the F layer in 
the FS bilayer calculated numerically for T„, = \/llT c . Plots (a)-(d) 
are calculated for h = kT c and (e)-(h) for h = 3nT c . The temperature 
T = 0.ir c . Parameters of the F/S interface are 7= Yb2 = 0.01 (a),(e): 
d f /% n =0.5; (b),(f): df/$ n = 1; (c),(g): d//|„ =2; (d),(h): = 3. 
The approximate analytical solution [5] is shown by dashed red lines. 

by proceeding to the analytical continuation in Eqs. (2), 
(3), (5), (6), (7) of the quasi-particle energy ico — > E + iO 
and using the A(x) dependence obtained in the previ- 
ous step, we find the Green functions by repeating the 
iterations until convergency is reached. We define the 
full DOS N{E) and the spin resolved DOS N t(l) (E), 
normalized to the DOS in the normal state, as 

N(E)=[N r (E)+N i (E)]/2, (8a) 
N m) (E) = Re [cos m) (ico -^E + i0)] , (8b) 

where 9/, is the boundary value of Of at x — —df/2. In 
case of a long F layer (df 3> it is also possible to 
obtain an analytical expression for the DOS at the free 
boundary of the ferromagnet [5,11]. 

In Fig. 2 and 3 we plot the DOS energy dependence 
for different df. Fig. 2 corresponds to the case of the 
absence of spin-flip scattering (infinite T,„) and Fig. 3 
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Fig. 4. (Color online) Current-voltage characteristics of a SIFS 
junction in the absence of spin-flip scattering for different values of 
the F-layer thickness df. The temperature T = 0.ir c . The exchange 
field h = (black line, which correspond to the case of a SINS 
junction), h = 7lT c (red dash-dotted line), and h = 3tcT c (blue dashed 
line), (a): d f /£, n = 0.5; (b): d f /^„ = 1, (c): d f /^„ = 2, and (d): 
df/i, = 3. 

to the case of x m = 1 / %T C . In both figures we plot the 
DOS for two chosen values of exchange field, h = nT c 
for plots (a)-(d) and h = 3nT c for plots (e)-(h). In Fig. 2 
we see that for any chosen h at small df the full DOS 
turns to zero inside a mini-gap, which vanishes with the 
increase of df. Then the DOS at the Fermi energy N(0) 
rapidly increases to the values larger than unity and with 
further increase of df it oscillates around unity while it's 
absolute value exponentially approaches unity. In the 
case of long enough ferromagnets we can observe DOS 
peaks at E — h [we notice that %T C s» 1.79A]. Also for 
df > 1.5 % n [5] we can use an analytical approximation 
for N(E) [5,1 1]. For smaller df it is incorrect to use this 
approximation. It is shown by the dashed red line and is 
in rather good agreement with the numerical solution. 
In the numerical curves the peaks at E = h are smeared 
because we used finite 7 = 0.01 for the transparent F/S 
interface at x = df/2. 

In Fig. 3 we observe similar tendencies. However for 
smaller x m the mini-gap closes at smaller df, the period 
of the DOS oscillations at the Fermi energy increases 
and the damped exponential decay occurs faster. Also 
the DOS peak at E — h is smeared. 

4. Current-voltage characteristics of SIFS 

We calculate the current- voltage characteristics 
(CVC) of a SIFS junction at T = 0.1T C using Eq. (1). 
Plots for different parameters are shown in Figs. 4 and 

5. For comparison we also present the CVC of a SINS 



Fig. 5. (Color online) Current- voltage characteristics of a SIFS 
junction for T„, = 1/kT c . The temperature T = Q.IT C . The exchange 
field h = (black line, which correspond to the case of a SINS 
junction), h = 7cT c (red dash-dotted line), and h = ?>nT c (blue dashed 
line), (a): d f /$, = 0.5; (b): d f /&, = 1, (c): d//$, = 2, and (d): 
d f l%, = 3. 

tunnel junction, i.e. a junction with a normal metal 
interlayer instead of a ferromagnet. We can see that 
for a certain range of parameters the CVC of a SIFS 
junction exhibits an interesting N-like feature [we can 
observe it for h = 3nT c in Fig. 4 (b) and for h = 7zT c in 
Fig. 4 (c)]. It corresponds to the case of a large subgap 
DOS in Fig. 2 (c),(f). For finite spin-flip scattering the 
features of CVC are smeared. At large df they totally 
disappear. A more detailed analysis of CVC of SIFS 
junctions will be presented elsewhere [12]. 

5. Summary 

To summarize, we calculated the quasi-particle DOS 
in the F layer of a SIFS junction in close vicinity of the 
tunnel barrier and use it to obtain the current-voltage 
characteristics of a SIFS junction. The developed for- 
malism may be used for estimations of dissipation in 
SIFS-based qubits, which will be done elsewhere. 
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